Specific heat and electronic states of superconducting boron-doped silicon carbide 
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The discoveries of superconductivity in the heavily-boron doped semiconductors diamonc 1 (C:B) in 2004 and 
silicon 2 ^ (Si:B) in 2006 have renewed the interest in the physics of the superconducting state of doped semicon- 
ductors. Recently, we discovered superconductivity in the closely related "mixed" system heavily boron-doped 
silcon carbide (SiC:B). 3 Interestingly, the latter compound is a type-I superconductor whereas the two afore- 
mentioned materials are type-II. In this paper we present an extensive analysis of our recent specific-heat study, 
as well as the band structure and expected Fermi surfaces. We observe an apparent quadratic temperature depen- 
dence of the electronic specific heat in the superconducting state. Possible reasons are a nodal gap structure or 
a residual density of states due to non-superconducting parts of the sample. The basic superconducting param- 
eters are estimated in a Ginzburg-Landau framework. We compare and discuss our results with those reported 
for C:B and Si:B. Finally, we comment on possible origins of the difference in the superconductivity of SiC:B 
compared to the two "parent" materials C:B and Si:B. 

PACS numbers: 74.25.Bt; 74.62.Dh; 74.70.-b; 74.70.Ad 



I. INTRODUCTION 

Diamond and silicon are wide-gapped semiconductors / in- 
sulators which exhibit indirect energy gaps of about 5.5 eV 
(diamond) and 1.1 eV (silicon). They are well-known for 
their outstanding physical properties and technical applica- 
tions, e. g. the excellent heat conductivity of diamond, its 
withstanding of high electric fields, or the numerous appli- 
cations of silicon in semiconductor technologies. It is well- 
known, too, that the physical properties of these and other 
semiconductors can be influenced by charge-carrier doping 
either by donor or acceptor atoms which changes their resis- 
tivity many orders of magnitude leading to intriguing prop- 
erties. Small doping concentrations are widely used in the 
application of semiconductors. At higher doping levels the 
systems undergo a semiconductor-to-metal transition above 
a certain critical doping level, i. e. charge-carrier concentra- 
tion, and further doping might even lead to superconductivity. 
From the theoretical and experimental point of view super- 
conductivity in doped semiconductors is an outstanding issue. 
The prediction of superconductivity in Ge and GeSi and the 
suggestion that other semiconductor-based compounds may 
also exhibit superconductivity at very low temperatures were 
given by Cohen already in 1964. 4 Indeed, some examples have 
been reported so far, e. g. self-doped Ge A TeJ^ Sn^ TeP do ped 
SrTiC^P or more recently doped silicon clathratesP^Eni The 
superconductivity of the doped silicon clathrates is the first 
example of compounds exhibiting superconductivity in a co- 
valent tetrahedral sp 3 network with bond lengths similar to 
those in diamond. 

However, before 2004 (C:B^ and 2006 (Si:B P supercon- 
ductivity was never reported for diamond and cubic silicon in 
the diamond structure although there are several studies avail- 



able concerning hole-doped induced metallicity in carbon and 
silicon using boron, nitrogen, or phosphorus, e. g. Refs. QT| 
[12] and the references therein. Therefore, it was an impor- 
tant progress to find superconductivity in these compounds 
upon boron doping, which attracted a lot of interest and stim- 
ulated ma ny theoretical and experiment al studies in the last 
four years! 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 2 ^ 23 ' 24 '^' Boron has one par- 
tially filled electron less than carbon or silicon and hence acts 
as an acceptor leading to hole doping. Both compounds are 
type-II superconductors with T c values of 11.4K (C:B) and 
0.35 K (Si:B). The upper critical fields are H C 2 ~ 8.7 T and 
0.4T,respectivelyP22l 

In order to explain the superconductivity in diamond 
theoretical studies point towards two different scenar- 
ios: (i) it is th e result of a simple electron-phonon 
interactior ] 14 ! 15 ! 17 ! 19 -! and (ii) it is caused by a resonating 
valence-bond mechanism ! 13 ! 27 ! 28 ! The former model is based 
on a conventional electron-phonon mechanism where the 
charge carriers are introduced into intrinsic diamond bands 
leading to a three-dimensional analog of the two-dimensional 
superconductor MgB2. The superconductivity is attributed to 
holes located at the top of the zone-centered cr-bonding va- 
lence bands which couple strongly to optical bond-stretching 
phonon modesP^l The latter model attributes the supercon- 
ductivity to holes in the impurity bands rather than in the 
intrinsic diamond bandsP^ With the premise that the dop- 
ing level in superconducting diamond is close to the Mott 
limit the randomly distributed boron atoms, i. e. their random 
Coulomb potential, may lift the degeneracy of the boron ac- 
ceptor states leading to a narrow half-filled band from which 
superconductivity develops. However, spectroscopical studies 
seem to support the former explanation and rule out the lat- 
ter suggestion] 29 ! 30 ! although a complete understanding of the 
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(a) 3C-SiC (b) 6H-SiC 




FIG. 1: (color online) (a) Unit cell of diamond-related cubic 3C-SiC. 
The three bilayers consisting of C and Si layers are emphasized. The 
stacking sequence is ABC - . . . The green arrow denotes the { 1 1 1) di- 
rection whereas the gray rods refer to the tetrahedral bond alignment 
of diamond. The cube defines the conventional unit cell of 3C-SiC, 
which consists of four formula units SiC. (b) Four unit cells of hexag- 
onal 6H-SiC. The cuboid defines one unit cell. The six bilayers of the 
stacking sequence ABCACB - . . . are emphasized. The green arrow 
denotes the (001) direction. One unit cell of 6H-SiC consists of six 
formula units SiC. For the drawings the software Vesta was usedPl 



superconducting phase is not yet obtained! 23 ! 28 * 31 ! Recently, a 
theoretical study suggested the possibility to achieve super- 
conducting transition temperatures on the order of 100K in 
C:B due to the exceptionally high Debye temperature of dia- 
mond and under the precondition that the doped boron atoms 
are ordered! 2 ^ 

In Ref. [3] we reported the discovery of superconductivity 
in a closely related system originating from a well-known and 
widely used semiconductor, namely boron-doped SiC, the sto- 
ichiometric "mixture" of the two afore discussed "parent" ma- 
terials. SiC is used increasingly for high-temperature, high- 
power, and high-frequency applications due to its high thermal 
conductivity, the existence of large band gaps, strong covalent 
bondings, chemical inertness, or its high tolerance to radiation 
and heat. Another hallmark of this system is the huge num- 
ber (abouiSl 200) of crystal modifications with cubic ("C"), 
hexagonal ("H"), or rombohedral ("R") symmetry of the unit 
cell! 33 l 34 lThey are usually referred to as wC-SiC, mH-SiC, and 
mR-SiC, respectively. The variable m gives the number of 
Si - C bilayers consisting of a C and a Si layer stacking in the 
unit cell. However, most of the available studies refer to the 
following polytypes: 34 cubic 3C- (zincblende structure, space 
group F43m (T?); "ordered" diamond) and hexagonal 2H-, 
4H-, and 6H-SiC (wurtzite, moissanite-4H, and -6H structure, 
all space group P63111C (Cf )). The 3C- (2H-) polytype is the 
only "pure" cubic (hexagonal) modification, all other mH-SiC 
polytypes consist of hexagonal and cubic bonds!^ The cubic 
3C-modification is also labeled as /J-SiC, whereas the hexag- 
onal polytypes are generally denoted as a-SiC. Fig. [T] gives 
a sketch of (a) the diamond-related modification 3C-SiC and 
(b) the hexagonal 6H-SiC. The C - Si bilayers are emphasized. 
In 3C-SiC both elements form face-centered cubic sublattices 
which are shifted by (1/4, 1/4, 1/4) with respect to each other. 



TABLE I: Basic parameters of 3C-SiC and 6H-SiC at room 
temperature The parameter V denotes the volume of the conven- 
tional unit cell, V mo \ = Vo ■ N^/t the molar volume where t is the 
number of formula units SiC in the unit cell ("f. u. / unit cell"). In 
our analysis we will use the average of SlC and V™f SlC because 
the sample used contains both polytypes; see text. 
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energy gap (eV) 


2.2 


3.02 


&D (K) 


1270 


1200 



Along the (111) direction the bilayer stacking in 3C-SiC is 
ABC -. . . For the polytypes 2H-, 4H-, and 6H-SiC it is along 
the (001) direction ABAB -. . . , ABAC -. . . , and ABCACB - 
. . . Some basic parameters of undoped 3C- and 6H-SiC at 
room temperature are summarized in Table[l] 

Depending on the crystal modification, pure SiC exhibits 
an indirect energy gap between ~ 2 eV (3C-SiC) and ~ 3.3 eV 
(2H-SiC)! 3 ^ Slightly doped SiC with donors and acceptors 
was intensely studied for nitrogen, phosphorus, boron, alu- 
minum, etc. by ion-implantation or thermo-diffusion doping. 
Compared with other dopants, boron was found to have a 
much faster diffusion rate in SiC. Diffusion processes medi- 
ated by the silicon interstitials and by carbon vacanci es have 
been proposed to explain such fast diffusion rates j2ZI2 g| 39 l 40 l 
Under silicon-rich conditions the carbon-site substitution is 
dominating. 41 Among other dopants, the insulator-to-metal 
transition was observed recently in nitrogen-doped 4H-SiC at 
carrier concentrations above lO^cm" 3 ! 4 ^ 

In this paper we report a specific-heat study on SiC:B 
and give a detailed analysis. Moreover, the density-of-states, 
band dispersions, and two- and three-dimensional plots of the 
Fermi surfaces for 3C-SiC are presented. We estimate the 
basic superconducting parameters and compare our findings 
with the reported results for C:B and Si:B. Finally, we com- 
ment on possible origins of the differences between the three 
superconducting systems. 
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H. EXPERIMENT 

The preparation and characterization of our samples is de- 
scribed in detail in Ref. [3j We studied several samples from 
different growth processes reproducing the general findings 
presented in this paper. The particular sample used in this 
study is identical to that used to map out the H-T phase di- 
agram in our previous study, namely "sample 1", referred to 
as "SiC-1" in this paper. The hole-doping charge-carrier con- 
centration of SiC-1 was estimated to be n — 1.91 • 10 21 cm' 3 
by a Hall-effect measurement.^ We note, that all of our so-far 
prepared samples are polycrystalline materials and not single 
phase. We found phase fractions of 3C-SiC, 6H-SiC, and Si. 
In spite of this result the residual resistivity for specimen SiC- 
1 is already as low as 60 jtOcm. The residual-resistivity ra- 
tio RRR = P3ook/Pi.5k amounts to 10. SiC-1 undergoes a 
sharp superconducting transition around 1 .45 K. The thermo- 
dynamic critical field is estimated to be ~ 115 Oe. In contrast 
to the type-II superconductivity in C:B and Si:B the nature 
of the superconductivity in SiC:B is type-I: We find a clear 
hysteresis in the temperature (field) dependence of the AC 
susceptibility between cooling (field-down sweep) and subse- 
quent warming (field-up sweep) runs indicating the hallmark 
of type-I superconductivity as discussed in Ref.|3] 

Specific-heat data was taken by a relaxation-time method 
using a commercial system (Quantum Design, PPMS). First, 
we applied a degaussing procedure before the measurement in 
order to reduce the remanent field of the magnet. Second, the 
addendum heat capacity was measured at Oe. Next, specific- 
heat data was taken in H — Oe and subsequently in 200 Oe, 
for which the addenda data was not measured, because the 
difference is expected to be negligibly small. However, this 
procedure lead to a small but visible artifact in the in-field 
normal-state data for 0.45 K < T < 0.6 K. Thus, the corre- 
sponding data points were removed and not used for the anal- 
ysis. 



III. RESULTS AND ANALYSIS 
A. Specific heat 

Fig. |2](a) summarizes the temperature dependence of the 
specific heat c p for specimen SiC-1 in the superconducting 
state (H = Oe) and in the normal-conducting state (achieved 
by applying a magnetic field H = 200 Oe > H c ). The raw data 
of this figure is the same as that used for Fig. 4 in Ref. [3]ex- 
cept the removed data points which were affected by an exper- 
imental artifact. SiC:B is a bulk superconductor as indicated 
by the clear jump of c p at T c . The inset of Fig. 151(a) shows 
the specific heat up to ~ 360 K for comparison. TTie room 
temperature value is still clearly below the classical high-T 
Dulong-Petit limit (49.88 J/molK). 

The solid curve in Fig. |2j(a) is a fit to the in-field data for 
0.6 K < T < 2 K applying the conventional Debye formula 

c P = c ph + c d = y„T + pT 3 (1) 




r(K) r(K) 



FIG. 2: (color online) Specific heat of SiC-1: The (red) • symbols 
denote the data in zero magnetic field. The (blue) a refer to data 
measured in a magnetic field H = 200 Oe > H c , representing the 
normal-state specific heat, (a) Specific heat c p as measured: The line 
is a fit to the in-field data using the standard Debye formula (Eq.[TJ. 
The inset shows the specific heat up to 360 K. (b) Electronic specific 
heat c c i/r: The lines are an entropy-conserving construction in order 
to estimate the intrinsic jump height; see text. 



with the Sommerfeld coefficient of the normal-state specific 
heat y n and the coefficient of the phononic contribution p 1 
as adjustable parameters. The fit yields a very good de- 
scription of the data below 2K. The obtained values are 
y n = 0.29mJ/molK 2 and £ = 0.02mJ/molK 4 . From the 
latter value we determined the Debye temperature 6>d using 
p 1 = (12/5)?tAWaA:b/@p with the number of atoms per for- 
mula unit N - 2, the Avogadro number Na, and Boltzmann's 
constant k^, yielding 6>d = 590 K. This is surprisingly low, 
only half of the value reported for undoped SiC (cf. Table|I|: 
0^ c * 1200K - 1300 K. We note that in the case of C:B a 
significant reduction of the Debye temperature to about 75 % 
of the pure diamond value (6>£ * 1860 K, @g ~ 625KjSlis 
reported, too. 45 In SiC:B this effect turns out to be even more 
pronounced. 

To further analyze the data we plot the electronic specific 
heat c e i - c p - c p h in Fig. |2j(b). The specific-heat jump 
starts slightly below T c * 1 .5~K coinciding with the results 
obtained by our AC susceptibility and resistivity measure- 
ments. However, the superconducting transition in c p is rather 
broad. The lines in Fig. |2j(b) are an entropy-conserving con- 
struction in order to estimate the intrinsic jump hight Ac e i. 
The "jump" temperature 7** = 1.31 K indicated by the per- 
pendicular (green) line in Fig. |2](b) is lower than the onset 
temperature T c reflecting the broadness of the transition. The 
jump height is estimated to Ac e i/y R T* * 1. The obtained 
value is only two thirds of the weak-coupling BCS expecta- 
tion, namely 1.43. For C:B and other semicondcutor-based 
superconductors, e.g. GenosTe, an even smaller jump height 
as low as 0.5 is reportedpSEH However, the overall shape of 
the C:B specific-heat data given in Fig. 3 of Ref. 45 is quali- 
tatively different compared to our data. The authors report a 
very broad transition consisting of two well-separated transi- 
tions. 

Next we focus on the question of the superconducting gap 
symmetry. We try the following two models to describe our 
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experimental data: 

Model (i) Assuming an isotropic gap structure: 
The simplest approach to obtain information about the su- 
perconducting gap is given by the conventional BCS text-book 
formula 47 c e \(T)/T cc exp(-A(0)/r)/7\ However, paying re- 
spect to the facts that on one hand the exponential behavior is 
only expected well below T c and on the other hand data below 
approximately 0.45 K is lacking leads to the idea to replace 
the exponential term by tabulated numerical specific-heat data 
calculated in the standard weak-coupling BCS framework,ES 
which is in principal valid up to T c . Therefore, we fitted the 
tabulated data with a polynomial c B , cs (15 th order) leading to 



ca.Cn/T cc cl cs (T)/T. 



(2) 



Next, considering that the samples are not single phase, it 
is reasonable to assume an additional T-linear term y res r 
reflecting a residual density of states originating from non- 
superconducting metallic inclusions. This modifies Eq.|2]as 
follows: 



cdT)/T = y res + y s ■ c» cs (r)/r. 



(3) 



Since the entropy related to a residual term y res T* does not 
contribute to the specific-heat jump, the prefactor of the BCS 
term is given by y s = y n - y res . Therefore y res is the only 
adjustable parameter in this approach. 

Model (ii) Assuming a power-law behavior of the electronic 
specific heat: 

At temperatures well below T c , a superconducting gap 
structure with nodes leads to the power-law behavior 



cdT)/T = y res + a ■ T h 



(4) 



with b — 1 or 2 for line or point nodesP^ES Here, we pay 
respect to a residual contribution, too. 

The results obtained by applying both models are summa- 
rized in Fig. [3] Applying Eq. ([3]) to the data yields the green 
dashed curve, Eq. Q the black dotted line. The estimated y 
factors for model (i) are given, too: y n (normal-conducting 
state), y s (superconducting state), y res (residual contribution). 

Model (i) Applying Eq. Q to the data for T < 0.7 K gives 
a reasonable description as seen in Fig.[3](green dashed line). 
The residual y coefficient amounts to y res = 0.14mJ/molK 2 . 
Above 0.7 K the fit undershoots the data corresponding to a 
distribution of T c values reflected in the broadness of the tran- 
sition. An entropy-conserving construction using the fit result 
instead of the linear approximation shown in Fig. |2](b) yields 
a slightly higher "jump" temperature T* = 1.35 K due to the 
downwards curvature around T c . Paying respect to the resid- 
ual contribution y les T* and evaluating the jump height with 
the Sommerfeld parameter of the superconducting part of the 
sample, y s = 0.16 mJ/molK 2 , gives almost the value predicted 
by the BCS theory: Ac el /y s T* = 1.48. 

Model (ii) The assumption of a power-law behavior in or- 
der to describe the data yields an even better description: As- 
suming a linear temperature dependence of c e \/T reproduces 
the experimental data in the whole temperature range below 
« 1 . 1 K, i. e. below the transition down to 0.45 K and extrapo- 
lates further down to K without any indication of a residual 
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FIG. 3: (color online) Electronic specific heat of SiC-1: The (red) • 
/ blue A symbols denote again zero-field / in-field data. The lines are 
fits to the data assuming an isotropic gap (dotted line) and a nodal gap 
(dahsed line). The dashed-dotted horizontal line denotes the value of 
the Sommerfeld parameter in the normal-conducting state y„. The 
two other y values are estimates related to model (i) and denote the 
superconducting y s and a possible residual Sommerfeld parameter 
y res ; see text. 



contribution y KS T* in contrast to the results obtained by apply- 
ing model (i) to the data. The resulting fitting curve is shown 
in Fig. [3] (dotted black line). It was obtained by adjusting only 
the prefactor a and fixing b = 1 and y res = in Eq. The 
attempt to include a residual term to the fit gave y res ~ 0. A 
T-linear behavior of c e \/T is expected in the case of a gap 
containing line nodes, but only well below T c where the su- 
perconducting gap is nearly independent of temperature. It is 
expected that at higher temperatures T — > T c the specific heat 
is affected by the reduction of the gap magnitude and therefore 
deviates from the linear extrapolation. 50 A complex balance of 
different effects is needed to cause an apparent linear temper- 
ature dependence up to T c . We note, that a T-linear behavior 
of c e \/T up to T c has been reported for e. g. the heavy-fermion 
compounds URu 2 Si 2 and UPt 3 pHI] 

The obtained jump height using a linear entropy -conserving 
construction, Ac e i/y n T* m 1 (Fig.[2](b)), is similar to the value 
predicte d theo retically for a superconductor with a nodal gap 
structure! 50 ! 52 ! We note that in such a model c e \/T should ex- 
hibit a rounded maximum at T c rather than the triangle-like 
peak used in the simplified entropy-conserving construction 
shown in Fig.|2](b). Therefore the jump height would be even 
smaller than 1 and T* slightly higher. 

However, in spite of the satisfying description of the data 
following model (ii) it is still necessary to obtain data down 
to several lOmK to clarify the true nature of the supercon- 
ducting gap. It would not be surprising if the specific heat 
of the multi-phase sample used consists of an additional T- 
linear term due to a residual y res as suggested by the result of 
model (i). Therefore the apparent power-law behavior of the 
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TABLE II: Normal-state and superconducting properties of SiC:B 
compared to those reported for C:B (Refs. [T] and l45l > and Si:B 
(Ref. |2j. Note that the highest T c (H C 2) for C:B reported so far is 
11.4K (8.7 • 10 4 Oe)P The asterisked "*. . . *" values are prelimi- 
nary because they depend on the value of po which we believe is still 
not the intrinsic one; see text. The coherence length £ in the case 
of the type-II superconductor C:B was estimated using Eq. j5} (i. e. 
from n, y n , and T c ) for better comparison with the type-I supercon- 
ductor SiC:B. The numbers given in parantheses are the published 
values from Refs. 2 (Si:B) and 45 (C:B) calculated with Eq. ([8} (i. e. 
from H c2 ). 





SiC:B 


C:B 


Si:B 


n (cm 3 ) 


1.91 • 10 21 


1.80 • 10 21 


2.80- 10 21 


y n (mJ/molK 2 ) 


0.294 


0.113 




£(mJ/molK 4 ) 


0.0193 


0.0007 




D (K) 


590 


1440 




*—**■' ei / / n - 1 q 


1 


0.50 




i-i,, ( 1 1 Orm^ 


*60* 


2500 
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RRR 


*10.0* 


0.9 


1.2 


r c (0) (K) (onset) 


1.45 


4.50 


0.35 


ff c (0) (Oe) 


115 






H x (0) (Oe) 


80 






H cl (0) (Oe) 


type-I 


4.2 • 10 4 


4000 


k F (nm -1 ) 


3.8 


3.8 




m* (m e i) 


1.2 


1.7 




T(0) (fs) 








v F (m/s) 


3.8 • 10 5 






f(nm) 




0.34 




?(0) (nm) 


360 


80 (9) 


(20) 


,1(0) (nm) 


130 


160 




*gl(0) 


0.35 


2(18) 





experimental electronic specific heat for 0.45 K < T < 1 . 1 K 
extrapolating to for T — » K is rather striking. 



B. Superconducting Parameters 

Together with the resistivity, Hall-effect, and AC suscep- 
tibility data published in Ref. 3 we are able to estimate the 
basic superconducting parameters. They are summarized 
in Table [IT] along with the derived normal-state parameters. 
For comparison the so-far known corresponding parameters 
for C:B and Si:B are listed, too. From the charge-carrier 
concentration 4 ' 5 (n = 1.91 ■ 10 21 cm~ 3 ) assuming a single 
spherical Fermi surface we obtain the Fermi-wave number 
kp = (3n 2 n) 1 ^ 3 = 3.8nm _I . The effective mass is evaluated 
as m* = Oh 2 y n )/(V mo ik^k F ) = \.2m e \ with the bare-electron 
mass m e \ and the molar volume 4 ^ V mo i- The Fermi velocity 
v F = Hk F /m* amounts to about 0. 1 % of the speed of light. The 
mean-free-path is estimated to be ( — Mp/(pone 2 ) = 14 nm 
with the elementary charge e. The superconducting penetra- 
tion depth amounts to /t(0) = y/m* /(/j.Qtie 2 ) - 130nm. The 
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FIG. 4: (color online) Ratio of H c /H sc according to Eq. IS): The solid 
line is a linear fit to the data; see text. 



coherence length is estimated using the BCS expressiorF-^ 

#0) = 0A8hv F /(k B T c ), (5) 



which yields £(0) 

m 



360 nm. Hence, the Ginzburg-Landau 
parameters k gl = 0.96^(0) /£(0) = 0.35 < 1/ V2, clearly 
placing SiC:B in the type-I regime. 

An independent approach to estimate the GL parameter kql 
is to start with the two phase lines given in the H-T phase di- 
agram (Fig. 5 in Ref.[3j, i. e. the supercooling field H sc and the 
critical field strength H c . Following the GL theory the analy- 
sis of the observed supercooling behavior in AC susceptibility 
provides an upper limit of the GL parameter. Upon decreasing 
an external applied magnetic field at constant temperature the 
superconducting nucleation field is given by H C 2- For a type-I 
superconductor this is smaller than the thermodynamic criti- 
cal field H c , which leads to the effect of supercooling if more- 
over the GL parameter satisfies kq L < 0.417. However, Saint- 
James and de Gennes showed that the nucleation of super- 
conducting parts of a sample sets in near to the surface at the 
surface nucleation field H c i = 1.695// c 2 = 1.695 V2/cgl#c, 
larger than H C 2 if the sample is placed in vacuum. 5455 - In real 
experiments, however, the onset of superconductivity will be 
observed at a field H sc , which is larger than the ideal super- 
cooling H c y. the experimentally derived supercooling phase 
line in an H-T phase diagram will satisfy the inequality 
H sc > H C 3. Hence, the "real" GL parameter has to be smaller 



than the "supercooling" k< 



J47l56l 



1 



Her 



1.695 V2 #c 



0.417- 



7^ 



(6) 



We would like to emphasize that the qualitative observation of 
supercooling in the field-dependence of our AC susceptibility 
data already confirms that the GL parameter has to be smaller 
than 0.417. 

We calculated the ratio of H sc and H c deduced from field- 
sweep measurements at several temperatures. The result is 
shown in Fig.|4] Following the procedure used by Feder and 
McLachlan, 56 i. e. extrapolating the data to T = T c , yields 
in our case kql < 0.3, which is even smaller than the afore 
reported estimate derived from «, T c , and y n , supporting the 
conclusion that SiC:B is a type-I superconductor. 
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FIG. 5: (color online) Calculated density of states (DOS) vs. energy 
for zincblende 3C-SiC. The inset gives an enlarged view of the en- 
ergy range near to the Fermi level Ep. The origin of energy is taken 
at the valence-band maximum without doping. The dotted lines in 
both panels mark the experimental value of the Sommerfeld parame- 
ter y n = 0.29 J/molK 2 = 0.5 /(eV unit cell) and the respective energy 
shift due to the boron doping AE = 0.56 eV assuming rigid bands; 
see text for details. 
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FIG. 6: (color online) Band-structure calculation of 3C-SiC. Spin- 
orbit coupling is included. The upper three valence bands are the 
heavy-hole (hh, green), the light-hole (lh, red), and the split-hole 
(sh, blue) band. Panel (a) summarizes the band-structure calcula- 
tion based on the DOS calculation given in Fig.[5] Panel (b) gives 
an enlarged view along the [110] direction from the Brillouin zone 
center T towards the K point. The dotted lines mark the energy shift 
AE = 0.56 eV estimated from the data shown in Fig.|5]and the corre- 
sponding Fermi-wave number & F = 0.5-fcp(r-K) = 0.5- y/TS/4-2n/a; 
see text for details. Panel (c) gives an enlarged view of the zone cen- 
ter at T. The bands are split due to the spin-orbit coupling. 



However, our results also imply that SiC:B is a dirty-limit 
superconductor because the coherence length is much larger 
than the mean-free path: £(0) » t. The Ginzburg-Landau pa- 
rameter for a dirty-limit superconductor 47 is given by atql = 
0.71 5/1(0) /£ which is » 1/ V2 in the sample used due to the 
small I value. The large £gl / small I is mainly caused by 
the residual resistivity po- Among our samples the residual re- 
sistivity varies from 60 ^iQcm to the order of mflcm, all of 
them exhibiting a type-I behavior in the AC susceptibility. 
With this experimental finding and keeping in mind that the 
so-far prepared crystals are polycrystalline multi-phase ma- 
terials, it is reasonable to assume that the intrinsic value of 
po (and hence €) could be much lower (larger) than even the 
60 uI2cm (14nm) found for the sample SiC-1. The quantities 
given in Table |H| which are related to the value of po are not 
very reliable and therefore asterisked "*. . . *". A decrease of 
the residual resistivity to a few ^tOcm would be sufficient to 
shift kq L below the critical value of 0.417 < 1/ V2 in accor- 
dance with our experimental finding of a supercooled type-I 
superconductor. 

C. Band structure of 3C-SiC 

Calculated density of states (DOS) and band structure data 
provide another possibility to determine an upper limit of the 
Fermi-wave number and hence the GL parameter kgl using 
the experimental value of the Sommerfeld coefficient and 




For simplicity, we will focus only on the 3C-modification 
of SiC. We approximate, that all three valence bands are free- 



electron like. Moreover, we assume rigid bands, i. e. the band 
structure is independent of charge-carrier doping. 

The electronic band structure of 3C-SiC is calculated 
within the local density approximation (LDA) to the den- 
sity functional theory. The all-electron full-potential linear- 
augmented-plane-wave method is used to solve one-electron 
Kohn-Sham equations. All the relativistic effects including 
spin-orbit coupling are included to every self-consistent-field 
iteration. The results for the total DOS and the band disper- 
sions are shown in Figs. [5] and [6] respectively. Two- and three- 
dimensional plots of the three Fermi surfaces corresponding 
to the upper three valence bands of 3C-SiC are given in Fig. [7] 
Therein panel (a), (b), and (c) display three-dimensional rep- 
resentations of the heavy-hole (hh), the light-hole (lh), and the 
split-hole (sh) bands. Panel (d) gives their cross sections.^ 

The DOS corresponding to the experimental y n value is in- 
dicated by a dotted line in Fig. [5] The inset gives an expanded 
view of the relevant energy range. Using Eq. d7), the molar 
volume y mo i and the volume of the unit cell Vb (cf. Table[I]) 
to convert the corresponding units, y n = 0.29 mJ/molK corre- 
sponds to DOS ~ 6 ■ 10 21 states /(eV-cm 3 ) = 0.5 states /(eV- 
unit cell). The energy shift of the Fermi energy due to the 
charge-carrier doping was estimated to AE = 0.56 e V also in- 
dicated by a dotted line in Fig. [5] 

The band structure of 3C-SiC is shown in Fig.|6](a). Panel 
(b) provides an enlarged view of the relevant bands near the 
F point. In panel (c) the effect of spin-orbit coupling at the F 
point is shown, which is about lOmeV. 

Next we plot the estimated energy shift AE into the band 
structure plot as marked by a dotted line in Figs.[6](a) and (b). 
The corresponding Fermi-wave number for the heavy-hole 
band was estimated to k F m 50 % of the distance from T to K in 
the fee Brillouin zone which equals VT8/4-27r/a. The param- 
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FIG. 7: (color online) Plots of the Fermi surfaces of 3C-SiC around 
the T point of the fee Brillouin zone (given in black) for a Fermi 
energy below the top of the valence band by -0.56 eV. Only three 
instead of six bands are shown due to the smallness of the spin-orbit 
splitting. Panels (a), (b), and (c) display three-dimensional represen- 
tations of the heavy-hole (hh), the light-hole (lh), and the split-hole 
(sh) band. Panel (d) contains their cross sections. The colors of the 
respective band plots are the same as used in Fig.[6] 

eter a denotes the lattice constant of 3C-SiC a = 4.3596 A (cf. 
Tablejljl. This yields an upper limit of kp < 7.6nirT 1 , which 
is more than double the value evaluated before (3.8 nirT 1 ) as- 
suming a single spherical Fermi surface. 

Possible origins of the discrepancy between the two esti- 
mates of kp could be the neglect of the 6H-SiC phase fraction. 
Moreover, the assumption of a single spherical Fermi surface 
is too simple as can be seen in Fig[7] It is not known how the 
heavy-boron doping modifies the real DOS and band struc- 
ture, either, leading back to the question if the superconduc- 
tivity in this material evolves from intrinsic or impurity bands 
as discussed in literature for C:B. 



IV. COMPARISON AND DISCUSSION 

After the analysis of the experimental data we now focus on 
the question why the "mixed" compound SiC is a type-I super- 
conductor whereas the "pure" parent compounds silicon and 
diamond exhibit type-II superconductivity upon boron dop- 
ing. Let us therefore briefly compare the obtained parameters 
of SiC:B with those reported for C:BPFor Si:B no specific- 
heat study is available so far, hence a comparison is not pos- 
sible for all parameters. 

The charge-carrier concentrations for all three specimen are 
comparable and on the order of 2 • 10 21 cm -3 . However, the 
temperature dependence and the absolute values of the resis- 
tivity are different: SiC:B turns out to be a much better con- 
ductor exhibiting a metallic p(T) for T c < T < 300 K with 



an RRR value of about 10, whereas the resistivity of C:B and 
Si:B decreases slightly above T c . For Si:B the slope of the 
resistivity becomes positive above ~ 50 K, for C :B a slightly 
positive slope is observed only above 200KpEHln both cases 
the resistivity is almost temperature independent resulting in 
RRR values of about 1 . 

The Sommerfeld parameter y n is somewhat smaller for C:B 
compared to SiC:B, the coefficient of the phononic contribu- 
tion is much smaller for C:B resulting in a higher Debye tem- 
perature in the latter case. This is not surprising since the De- 
bye temperature of pure diamond is much higher than that of 
pure SiC. For both compounds the jump height of the specific 
heat (Fig.[2](b)) is much smaller than the BCS expectation for 
a weak-coupling superconductor. The superconducting pene- 
tration depths are similar for C:B and SiC:B (A(Q) * 150 nm) 
but the coherence lengths make the essential difference. The 
published values from Refs. 2 (Si:B) and 45 (C:B) have been 
estimated from the upper critical field strength H C 2 using the 
GL expression 

£ = V®o/2>rff c2 (0) (8) 

with the flux quanta (Do. Applying this formula the coherence 
lengths of C:B and Si:B are both on the order of 10 nm. These 
values of £ are given in parantheses in Table [III ^ n me case °f 
C:B we calculated £(0) using Eq. |5]) for a better comparabil- 
ity with SiC:B, too. The latter yields £(0) = 80 nm, whereas 
for SiC:B £(0) amounts to 360 nm resulting in different GL 
parameters kq L = 0.35 for SiC:B and 2 for C:B (which is 18 
using the published value of £ = 9 nm), cf. Table [II] Hence, 
SiC:B is a type-I and C:B a type-II superconductor. 

At the current state of research we can only speculate about 
the physical reasons for this different nature of superconduc- 
tivity in C:B / Si:B and SiC:B. In the case of C:B one appar- 
ent reason leading to a smaller coherence length and hence a 
larger GL parameter is the higher critical temperature of this 
superconductor: £(0) oc J*" 1 . However, this argument does not 
hold for Si:B, T c of which is much smaller than that of SiC:B. 
One can argue that SiC:B is a much cleaner system than C:B 
and Si:C. Hence, the coherence length ^(0) of Si:B (thin film 
and diffuse doping) might be limited by a very short mean- 
free path £ and thus the GL parameter is larger, leading to the 
speculation that "clean" Si:B could be a type-I superconduc- 
tor, too. 

Finally, we would like to mention a couple of apparent dif- 
ferences between the systems: 

Si-C bilayers: SiC is in a certain sense a "layered" sys- 
tem consisting of Si-C bilayers. Many polytypes are known 
distinguished by the stacking sequence of these bilayers in the 
crystal unit cell. In this sense one may refer to SiC as an "or- 
dered system". Our results 3 suggest that boron is introduced 
only into the carbon sites in SiC:B and hence only half of the 
crystal sites are directly affected by the disorder due to the 
hole-doping process, whereas in C:B and Si:B in principle all 
sites can be randomly involved. 

structure: For silicon and diamond the cubic crystal struc- 
ture seems to be important or a precondition for the appear- 
ance of superconductivity upon doping. In SiC the situation 
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is different. The multi-phase crystal used in this study con- 
tains cubic 3C-SiC and hexagonal 6H-SiC. At the moment 
we cannot rule out the possibility that both phase fractions 
contribute to the superconductivity which would be a clear 
difference compared to the two parent compounds. Moreover, 
Cohen suggested in 1964 that hexagonal SiC could exhibit 
superconductivity. 58 However, the same author predicts that 
most of the many-valley semiconductor based superconduc- 
tors should be type-II rather than type-IpJ 

band structure: In contrast to cubic diamond and sil- 
icon the "mixed" compound SiC breaks inversion symme- 
try. In crystals with inversion-symmetry, states with different 
spin orientations are degenerated. This is not true in general 
for crystals with broken inversion symmetry. The degener- 
acy might be lifted by the spin-orbit interaction. In 3C-SiC 
(zincblende structure) the degeneracy is preserved only along 
the [100] direction. Along e.g. the [110] direction (i.e. from 
the F point to the K poi nt in the Brillouin zone) the states with 
different spins split up! 59 l 6 ° l Using the above estimate of the 
Fermi-wave number for the heavy-hole band kp < 7.6nrrT 1 
we can give a rough estimate of the spin-orbit splitting for this 
band along F-K in 3C-SiCP A so = 0.02 me V, which is a 
rather small value. We note, that the light-hole, split-hole, and 
(in the case of electron doping) the lowest-conduction band 
exhibit larger values. 60 Using as an example kp = 3.8nm _1 
gives the same spin splitting because Aso of the heavy-hole 
band is almost constant in the interval 0.25 ■ 2n/a to 0.85 • 2n/a 
along the [110] direction of the Brillouin zoneP^ 

inversion symmetry: Moreover, a broken inversion sym- 
metry is known to give rise to a highly interesting nature of 
the superconducting ground state, including a spin singlet - 
triplet mixture! 61 * 62 ! In SiC:B we do not expect any uncon- 
ventional scenario based on the broken inversion symmetry 
because of the comparably light elements silicon and carbon 
without strong electron-electron interaction. 

charge-carrier concentration: In diamond, cubic silicon, 
and 3C-SiC the indirect band gaps are between the zone center 
(r point) and the X point of the Brillouin zone. For all other 
SiC polytypes the valence-band maximum is located at the F 
point, too, but the conduction-band minimum differs. For 6H- 
SiC it occurs at the M point! 36 ! 63 ! 64 ! If one takes into account 
spin-orbit interaction, the splitting of the band structure of di- 
amond and silicon is not affected, but for 3C-SiC the splitting 
will cause a shift of the valence-band maximum. 65 However, 
the boron doping in SiC removes electrons from the valence 
bands and therefore the difference in the semiconducting gaps 
might be of minor relevance. 

Nevertheless, it underlines again the importance of answer- 
ing the question whether the holes induced by boron doping 
in SiC reside in the intrinsic bands or form an impurity band, 
i. e. what is the nature of the metallic ground state, from which 



superconductivity develops? 



V. SUMMARY 

In summary we present a specific-heat study of heavily 
boron-doped silicon carbide SiC:B using the same crystal 
used in our recent publication reporting the discovery of su- 
perconductivity. In contrast to the type-II superconductivity 
in the two parent compounds, boron-doped diamond C:B and 
boron-doped silicon Si:B, the bulk superconductivity in SiC:B 
is type I. This is reflected in rather different values of the su- 
perconducting coherence length, i.e. 360 nm for SiC:B and 
only 80 nm for C:B, whereas the penetration depths are of 
the same order of magnitude. We presented two different ap- 
proaches to describe the data: assuming (i) an isotropic gap 
structure and (ii) a power-law behavior. The electronic spe- 
cific heat in the superconducting state is well reproduced by 
the former assumption with a residual density of states or by 
the latter assumption of a quadratic temperature dependence. 
The specific-heat jump height Ac e \/T at T c is about 1 or even 
smaller in the latter model and hence far away from the expec- 
tation in a BCS framework. However, due to the lack of data 
points below 0.45 K it is difficult to give a final conclusion 
about the superconducting gap structure. To further clarify 
the gap structure a specific-heat study in a dilution refrigera- 
tor system is desired. 

The origin of the different nature between the type-II su- 
perconductors C:B and Si:B on the one hand and the type-I 
superconductor SiC:B on the other hand remains at this state 
of research unclear. To clarify this intriguing issue further 
experimental and theoretical work is needed. From the ex- 
perimental point of view single crystalline samples are highly 
desirable. Moreover samples with only one phase fraction, 
either 3C-SiC or 6H-SiC, are eligible to answer the question 
which phase fraction is liable for the occurrence of supercon- 
ductivity in SiC:B. This work is currently under way. 
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